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We consider the problem of whether the canonical and microcanonical ensembles are 
locally equivalent for short-ranged quantum Hamiltonians of N spins arranged on a d- 
dimensional lattices. For any temperature for which the system has a finite correlation 
length, we prove that the canonical and microcanonical state are approximately equal on 
regions containing up to 0{N^^ (d+i)) spins. The proof rests on a variant of the Berry-Esseen 
theorem for quantum lattice systems and ideas from quantum information theory. 


I. INTRODUCTION 


In statistical mechanics there are two main ensembles (at zero chemical potential) that can 
be used to compute equilibrium properties of large systems: the microcanonical and canonical 
ensembles. Roughly, the first describes the physics of a system that is isolated and has total fixed 
energy The second describes the physics of a system that is at thermal equilibrium with a large 
environment at fixed temperature. It turns out that in many cases, although not all, the two 
ensembles give the same predictions for very large systems. There is a long sequence of studies 
aiming at elucidating under what conditions the two ensembles can be used interchangeably (see 
e.g. ||l}^6l and the discussion below). 

In textbooks the canonical ensemble is commonly introduced by considering the microcanon¬ 
ical ensemble of the system and a large environment and restricting to observables acting on the 
system only Under the assumption that the interactions of the system and environment are very 
weak, the canonical ensemble can be derived. However in many situations the assumption of 
weak coupling is not justified. For example, in many closed quantum systems small regions ther- 
malize (Zllll; in this case the remaining of the system is acting as an environment. It is therefore 
an interesting problem to find more general conditions that guarantee the equivalence of the two 
ensembles. Our main goal is to give one such condition: We show that short ranged interactions 
and a finite correlation length lead to the equivalence of ensembles for every sufficiently large fi¬ 
nite volume. The condition of a finite correlation length (and short ranged interactions) is known 
to be required (see e.g. |91). 


II. RESULTS 


We let A = {1,... ,n}'^ be a finite collection of vertices or lattice sites in d dimensions with 
= |A| = sites. We consider local Hamiltonians, acting on the Hilbert space PL = 
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dim T-Li = D, given by 


H = = (1) 

ieA U 

where we assume that the Hi are bounded, ||i^i|| < 1, and local in the sense that Hi acts only on 
sites j with dist(z, j) < k (for the Manhattan metric dist(.,.) in the lattice). 

For such A:-local Hamiltonians, we let pT '■= jZiT) be the canonical state at temperature 

T (also known as Gibbs state or thermal state) and Zi^T) ;= tr(e“^/^) the partition function 
(we set Boltzmarm's constant to unit). In the canonical ensemble at temperature T, averages are 
computed using pT- The energy density at temperature T is given by 

u{T) := ^tr(TFpT), (2) 

the specific heat capacity at temperature T by 

^ (tr[^VT] - {ir[HpT]f) , (3) 

and the entropy density by^ 

s{T) := -^tr[pTln(pT)]. (4) 


c{T) : = 


du{T') 


dT' 


Given regions X,Y C A, we denote by tr^^^^j^ the partial trace over the complement of X in A 
and for states p G V{H) (the set of density matrices acting on H), we denote pxY = trA\(v:uv)(/’)- 
Given two states p, a, their trace-norm distance is 


Up — (j||i := tr(|p — cr|) = max 2tr(M(p —u)) (5) 

0<M<I 


and quantifies how distinguishable the two states are. 

We say a state p G V{H) has (^, z)-exponentially decaying correlations (or a (^, 2 ;)-finite cor¬ 
relation length) if there are > 0 and z > 0 such that for every two regions X,Y C A with 
dist(X, Y) > 0, 


corp(X,y) 


max 


P,Q 

supp(P)cX 

supp((5)cv 


|tr((P (g) Q){pxY - Px® Py))\ 

ll^’IIIIQII 


< pjz^-disi{X,Y)/^^ 


( 6 ) 


where 


dist(X, y) := min dist(x,y). 

x£X,y£Y 


( 7 ) 


Given e G R and (5 > 0, let 


Me,5 := {ly : \E^ - eN\ < SVn}, 


^ Throughout, we denote by In (log) the logarithm to the base e (2). 


( 8 ) 
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and define the microcanonical state of mean energy e and energy spread (5\/iV by 

1 


^e,S • — 


I Me,51 




( 9 ) 


In the microcanonical ensemble averages are computed using Te,^. 

Finally, for a Gibbs state corresponding to a Hamiltonian as in Eq. 0 and with (^, z)- 
exponentially decaying correlations, we define (see Lemma |^: 


^k,i,z,T '■= Cd 


;max{fc,g}(z + 1)) 

./rW) 


2d 


max 


max{A;, + 1) In(iV) ’ T^c{T) 


( 10 ) 


where > 1 is a constant which only depends on the dimension of the lattice A. 

We can now state the main result. It shows that for general quantum many-body systems at 
non-critical temperatures (meaning that the canonical state has a finite correlation length), the 
canonical ensemble gives essentially the same predictions as the microcanonical ensemble, for 
every observable which acts on sufficiently small regions. 


Theorem 1. Let Ci be the set of all hypercubes contained in A = {1,... ,n}^ with edge length I G N, 
1 < I < and let N = nf > 2. Let the canonical state pT (corresponding to a k-local Hamiltonian 
as in Eq. 0) with energy density u{T) and specific heat capacity c{T) have (^, z)-exponentially decaying 
correlations. Let the microcanonical state mean energy e such that 


|e-u(r)| < ^fc{T)T‘^/N 

and energy spread 6y/N such that 

28Afc,^,,,rv'c(r)r2—^ < 5 < ^/c(T)TA 

Let e > 0. If 

56y^Afc,5,^,Tln2'^(iV) + (5 + e^)ln(iV) 2^\n{D)l^ + I + 2 ^ / eN 

eln(2) + ^ln(2) “ Un''(4)e^y 

then 

E \\{Te,s)c - iPT)c\\l < 7^/e, 


( 11 ) 

( 12 ) 

(13) 

(14) 


where the expectation is taken uniformly over Ci. 


We note the following: 


1. Eq. (131 is satisfied whenever N is sufficiently large and 


2. We do not need to take the average over regions C G if we assume the Hamiltonian is 
translation invariant. 

3. The condition of a finite correlation length is necessary. Indeed the two ensembles differ 
in the Ising model approaching the critical point, when the correlation length diverges, for 
regions of size 0{log{N)) (see e.g. IfJ). It is an open question if a similar result can be 
obtained for critical systems and small enough regions, assuming that correlations decay 
algebraically; in our proof it is important that the correlations decay exponentially as in 
Eq. 
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4. Any system is expected to have a finite correlation length whenever it is away from a crit¬ 
ical point. One can rigorously show that one-dimensional systems always have a finite 
correlation length at any temperature ITOl . while in any dimension there is a critical tem¬ 
perature (depending only on the geometry of the lattice) above which every system has a 
finite correlation length lUTI . 


An important step in the proof of the theorem will be to establish the following proposition, 
which we believe is of independent interest and which we prove in a stronger version (Propo¬ 
sition 1^ in Section V A| It shows that two states r and p are locally equivalent whenever their 
quantum relative entropy 


'S'(rllp) = tr(r(logr - logp)) (15) 

is 0(A^/('^+^)) and p has finite correlation length. 

Proposition 2. Let Ci be the set of all hypercubes contained in A = {1, ..., with edge length I G N, 
1 < I < and let N = n‘^ > 1. Let e > 0, the states p, r, and I such that p has (^, z)-exponentially 
decaying correlations and such that 

2{ln(P)i- + i + 2 rii 

e ^ fln(2) '-Vlii''(4xd ' 

Then 

E lire - pell 1 < 7Ve, (17) 

CeCi 

where the expectation is taken uniformly over Ci. 

We note that the state p does not need to be a thermal state, we merely demand it to have (^, z)- 
exponentially decaying correlations. If we do assume it is a thermal state then, as TS{t\\pt) = 
Ft{t) — Ft(pt), where, given a Hamiltonian H and temperature T, the free energy of a state r is 
given by 


FT(r)=tr(ifr)-r5(T), (18) 

the proposition shows that for temperatures away from criticality, any state of small free energy 
must have approximately thermal averages for local observables. Theoremj^will then follow from 
Proposition and showing that has small free energy whenever px has a finite correlation 
length. 


A. Beyond Microcanonical States 

How crucial is the use of the microcanonical ensemble? Proposition shows that not only the 
microcanonical state, but any state of small enough free energy is locally thermal. It turns out that 
Theorem[^can be extended in two additional respects: (1) It applies to any state that lives in the 
microcanonical subspace and has sufficiently large entropy. (2) Following [121 [l^, it applies to a 
generic state in the microcanonical subspace, with overwhelming probability with respect to the 
Haar measure. 
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Corollary 3. Let Ci be the set of all hypercubes contained in A = {1,. .., n}^ with edge length I G N, 
1 < I < and let N = Tif > 2. Let the canonical state pt (corresponding to a k-local Hamiltonian 
as above), with energy density u(T) and specific heat capacity c{T), have (^, z)-exponentially decaying 
correlations. Let ^dh e, 5 such that Eqs. 0 and ^12^ hold true. Let e > 0 and 

ln2'^(iV) + (5 + ez) ln(N) 2^ ln{D)l<^ + 1 + 2^1 
eln(2) ^ln(2) “ 2 


eN 


ln"*(4)^« 


1 

d+1 


(19) 


1. Let T a state on the subspace spanned by {\u)},^^Me s ^dh entropy 

5{r)>log(|M.,|)-i(^)”. (20) 

Then EceC; Ike - (pt)c||i < 7-v/e. 

2. Let T be a pure state drawn from the Haar measure on span{\u)}y(zM,,s. Then, with probability at 
least 1 — 


with 


IE Ike - ()Ot)c 111 < 7\/e + ?7 + 

CeCi 


pd 


n3/2 


p := 18^/^7rexp 



2yc(r) + 2 \ 

7n ) 


( 21 ) 


( 22 ) 


Here, the expectation is taken uniformly over Ci. 

The second part of the corollary is a direct consequence of the first part, the quantum Berry- 
Esseen bound in Lemmaand the result of lll^ll^ . which shows that a generic state in a energy 
subspace has the same local reductions as the microcanonical state. 


III. COMPARISON WITH PREVIOUS WORK 


The problem of equivalence of ensembles has been considered since the foundational work of 
Boltzmann and Gibbs. See I1T4I for a historical perspective. An intuitive explanation for the equiv¬ 
alence at non-critical temperatures is the following: Whenever there is a finite correlation length, 
the heat capacity, which determines the energy variance density, is of order 0(N), implying that 
the standard deviation of energy is of order 0{Vn). As energy is an extensive quantity (i.e. of 
order 0{N)), the distribution of energies in the Gibbs state is highly concentrated around its mean 
value for a large number of spins, and so one might be tempted to believe that only the micro- 
canonical subspace is relevant. However it turns out that this argument is too simplistic. Indeed 
it is easy to see that for any <5 > 0, Te^s arid px (with e = u{T)) are nearly orthogonal for suffi¬ 
ciently large N. Therefore any meaningful argument for the equivalence of ensembles must go 
beyond the distribution of energies and in some way restrict the kind of observables considered 
(for example, considering observables acting in small regions). 

The most fruitful direction explored so far has been to consider systems in the thermodynam¬ 
ical limit. In this regime one can prove the equivalence of ensembles on the level of thermody¬ 
namical functions ||TH3]-(5l (showing that the thermodynamical limits of the entropy density in the 
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microcanonical ensemble is the Legendre transform of the limit of the free energy density). One 
can also show it on the level of states, as we do here, both for classical |2| and, only very recently, 
for quantum systems |6i|. However the price of considering the thermodynamical limit—instead 
of the physically relevant regime of very large but finite sizes—is that no finite bounds can be 
obtained on the size of the regions on which the canonical and microcanonical states are close. 

In this respect Theorem]^ goes beyond the earlier work in several aspects: 

• It covers the general case of non translation-invariant models. 

• It is based on the assumption of a finite correlation length, which is simpler and more phys¬ 
ical than the assumption of a unique phase region employed in lH-El. 

• It gives explicit finite size bounds; for quite big regions of order the two en¬ 

sembles already look the same. 

• It shows that the equivalence holds true even for microcanonical states with very small 
energy spread, of order 0(log^'^(iV)) and substantially smaller than the value 0{y/N) that 
could have been expected. 

• It covers more general states than the microcanonical, showing that the important condi¬ 
tions are that the state has small free energy or that it is concentrated around a fixed energy 
and has sufficiently large entropy. 

• It shows that that any two microcanonical states arid Te'^s' are locally equivalent when¬ 
ever |e — e'\N < 0{'/N) and 0{log^‘^{N)/y/N) < 6 + 6' < 0(1) (assuming pT has a finite 
correlation length). 

It is an interesting open question to determine how small 5 can be taken. We note that the 
eigenstate thermalization hypothesis (ETH) states that even for <5 = 0, i.e. for a single eigenstate, 
one should already have the same local expectations values as the corresponding microcanonical 
state |[T5| . However, while believed to hold true for several systems, there are known counterex¬ 
amples to ETH, e.g., systems with many-body localization. 


IV. PROOF OUTLINE 


Our proof can be seen as a finite-size version of previous results HHSl relating the (mi- 
cro)canonical ensembles in the thermodynamical limit, in particular the recent work of Muller, 
Adlam, Masanes, and Wiebe |6J, who showed the equivalence on the level of states for quantum 
systems in the thermodynamical limit. There the authors obtained the result from two observa¬ 
tions, which we now briefly explain. 

Given a sequence of of translation-invariant Hamiltonians acting on finite volumes 
(with n spins) with a well-defined thermodynamical limit, we define the (Hermholtz) free energy 
density as 


/(T) := inf{/r(a;) : cj translation-invariant state}, (23) 

with /t(w) = u{uj)—Ts{u) the free energy density for a translation-invariant state ui in the infinite 
lattice limit. In 0 one is interested in the so-called one phase region, in which there is only one state 
(jj achieving the minimum in Eq. (j^, given by the KMS state associated to the sequence of finite 
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volume Gibbs states. This uniqueness condition holds true if the finite volume Gibbs states have 
a finite correlation length. The first observation of the authors of |6i| is that 


lim —Ft 

n—>oo n 


u{T),o{VN) 


lim — 

n^oo Ti 



(24) 


i.e., the free-energy density of the microcanonical ensemble converges to the free-energy density 
of the canonical ensemble (this fact is attributed to HU). The second observation is that because 
of the uniqueness assumption in Eq. ( |^ , '^u{T) o{Vn) converge to the same state and 

therefore for any fixed region A: 


lim 

n^oo 




hA„\A(PT") 


= 0 . 


(25) 


The proof of Theoremj^will have a similar structure to the argument above. Indeed Proposi¬ 
tion 1^ shows that every state of small enough free energy has its local reduced density matrices 
equal to the ones of the canonical state (if the latter has a finite correlation length). This is a finite- 
size analogue of the uniqueness of the minimizer in Eq. (231 (which, as we mentioned, can also 
be derived from the assumption of a finite correlation length). To prove Theorem]^ we show in 
Lemma 13 that 


S{TuiT),s\\pT) < Oilog^^N)). 


(26) 


This is a finite-size analogue of Eq. ( p4) | and follows from a version of the Berry-Esseen Theorem 
for quantum lattice systems that we prove in Ref liTTl (see Lemma|^. Note that since the maximum 
value of S{t\\pt) is 0{n) (as the reference state is a thermal state), the Eq. already suggests 
that the two states are not very different. Theoremj^then follows directly from Propositionand 
Eq. @. 


We now give a quick summary of the argument behind the proof of Proposition The proof 
in earnest is given in Section V A We use four variants of the quantum relative entropy. The first 
is the quantum Kullback-Leibler divergence defined before by: 


5(r||p) = tr(r(logr - logp)). 


(27) 


We also use its smoothed version: 


5*'(t||p) := min S’(f||p), (28) 

^ 6 Ss(r) 

with Bs{t) := {f : ||t — f ||i < e} the set of states that are e-close to r. We also consider the 
max-relative entropy of two states r and p KlSl : 

*S'max(T||/o) := {minA : r < 2Ap}, (29) 

and its smooth version 

‘S'max('^llp) := min S'max('rllp). (30) 

The relative entropies are related as follows 




< S'^{t\\p) < 


5(r||p) -h 1 


+ log 


e 


1 — e 


(31) 
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The second inequality is known as quantum substate theorem |[T9l|20l. 

For simplicity in this proof sketch we consider the one-dimensional case d = 1, N = n, 
leaving the general case to the actual proof. The set of all intervals of length I is given by 
Cl = {Cl,Cn-i+i} with Ci = + + l -1}. Thus 

E \\tc-Pc\\i = -r—r Y] \\rCi - PCi\\i- (32) 

ceCi n-l + 1 ^ 

2=1 

We may now group the sets Q such that the Q within each group are separated from each other 
by a distance of r. There are at most I -|- r such groups and within each group are m sets Ci. 
Let us now focus on one group and let Ci,..., Cm the sets in this group. Eqs. ( [T^ , plj l, and the 
monotonicity under partial trace give 

S^irci...Cm\\PCi...Cm) ^ Yn- (33) 

Lemma shows that if p has finite correlation length ^ then 

\\pCi...c,n - PCi®---® PCJ\i< (34) 


which can be made arbitrarily small by increasing r. By the data processing inequality for the 
smooth max-relative entropy HTSl (see Lemma in Section VD|| and Eq. ^3\ one then has 


5'max('rCi...C™l|PCi ®---®PCm)^ Yn- 


(35) 


with e' = 2y/e + VSmn^^. From Eq. (j^, we find that the above bound also holds for 
5'^ {TCi...Cm\\PCi ® ■ ■ ■ ® PCm)- Then by subadditivity of the von Neumann entropy it follows that 


^ Y{TCi IIpcJ < v^- (36) 

2=1 


Pinsker's inequality gives 


\Js^'{TCiWPC,) > \\TCi - pCih - e' 


(37) 


and hence, due to our choice of I < -^/en given by Eq. (16 ' and since by construction m n/{I+ r), 
we have 


E lire - pcWi <£' + (en)^/^ < Y + VY j. 

C&Ci n-l + l ^ 

and the result then follows from setting r ~ 


e/n, 


(38) 


V. PROOFS 


A. Proof of Proposition]^ 

Let 1 < I < n and Ci the set of all cubes in A with edge length I, i.e.. 

Cl = ^Ci I i G {1, n — I j- , Ci = i {0 ,... ,l — l}'^. 


( 39 ) 
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l + r 


1 2 ■■■ I ■■■ 


i 


n — I 1 
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00 
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00 
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00 


OOOOOOOOOOOO 

000000 

000000 

l + r 

00000 o_o 00000 

0 0 o_o 0 0 
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00 
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0000 
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00 
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0 0 
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00 
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t t t t 


*1 n + / + r • • • *i + (^ + r)ji • • • ii + (I + r)(m - 1) 


*2 


*2 + (i + »')(?n - 1) 


FIG. 1 : The lattice A = { 1 ,..., n}‘^ for d = 2 . The set of all cubic subsets of edge length / is C; = {Ci C A | / G A;}, 
where Ct = i + { 0 , 1 }'^ and A; = {l,n — I + l}'^, which we decompose as A; = lJiG{i i+r}'^ ■ Blue 

sites are the elements of for / = 4 , r = 2 , and i = (ii,J2) = ( 2 , 3 ). Blue squares are the corresponding 
cubic subsets of edge length I, Ci+(^i+r)j, which are separated by r sites. Black sites are the elements of for 
i = (11,12) = ( 1 , 1 ) and black squares indicate the corresponding Here, j € { 0 ,... ,to — l}'^ and 


m = \ 


n—Z+1 
Z+r 


1- 


Writing A; 


{l,n 


/ + 1 }'^, we thus have 


E \\tc-Pc\\i = 

C&Ci 


(n 


1 

l + lY 


ikci - pciWi- 

ieAi 


( 40 ) 


We now decompose A; into a union of sets such that for every j, j' G j + f, we have 
dist(C'j, Cf) > r as in Fig.j^ Define m := and 


such that 


and 


= u {* + (^ + r)j} n A; 

f S{0,...,m— 

^ |, + (; + r)j e A; I j e {0,..., m - l}*^! 
A,= (J £>''>, |A,|= ^ 


E \\tc - pcW 

C&Ci 


1 


E E 


\\TCj - PCjWl- 


( 41 ) 


( 42 ) 

( 43 ) 


We now need the following lemma. 

Lemma 4. Let 0 < e < 1 and Ci,..., Cm C A with dist{Ci, Cj) > Ofor i / j. Let r, p states such that 
p has {Y z)-exponentially decaying correlations and let 

M 

i =2 
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Then there is a state -rcx-Cm 

5'max(vrCi...CMll/5C’i ® ' ' ' <8) pcj < ^ + log (f~) + (i^) 

and \\ttcvCm “ '^Ci-CmIIi < 2^ + 

Proof. The quantum substate theorem 11191 in the version of Ref. 120 1 implies the existence of a 
state ^Ci-Cm such that 

W^Ci-Cm -'TCi-CmIIi ^ 2^ (46) 


and 


5'max(7rCi...CMll/5Ci-CA/) < 


S{tCi-Cm\\PCvCm) + 1 


+ log 


1 


1 - e 


such that, using the monotonicity under partial trace of the quantum relative entropy, 

1 


5(r||p) + 1 

^max{TrCi-CM\\PCi-CM) ^-^ log 


1 - e 


=: A. 


This proves the statement for M = 1. For M > 1, we use Lemma to find 

M 

IIpCi ®---®PCm - PCi-CmWi < =: c. 

i=2 

By Lemma]^ if k := 2^c < 1 then there is a state vrci-.-CAf such that 

5'max(vrCi-CM WPCi ® ® PCm) < ^ + log f 


(47) 


(48) 


(49) 


(50) 


and \\ttci-Cm ~ ^Ci - CmIIi — which, in combination with Eqs. (46' and (50 ', proves the 

statement for M > 1 as by the triangle inequality \\t^Ci-Cm ~ ^Ci-Cm IIi — 2-\/e + v8k. □ 

We are now in the position to prove the following stronger version of Proposition]^ 

Proposition 5. Let e > 0, the states p, r, and / G N, 1 < Z < such that p has z)-exponentially 
decaying correlations and such that 


^((2" - 1) ^ 2^ 


5(r||p) + 2 


< e(n-Z + l)“ 


(51) 


Then 

T. \\tc- Pc\\i<(V2 + 2 + ^Ji^)V2e. (52) 

CeC; V / 

Here, [•] denotes the smallest integer not less than ■ and W the solution to z = W{z)e^^^\ W{z) > 0 
(one of the real branches of the Lambert W function). 

Proof. We set out to combine Lemma with the following basic properties of the quantum relative 
entropy and trace norm. 
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[a] (Pinsker's inequality) \\p — (j\\\ < ln(4)S'(/9||cr), 

[b] (Relation with 5max US) S{p\\a) < Sm^^{p\\a), 

[c] (Super-additivity^) SiT^A,\\pAj) < S{7rA^-AM\\pAi ® • • • 0 pam)> 

[d] (Monotonicity under partial trace) ||pa — u^Hi < ||p — cr||i for all A c A. 


To this end, let 0 < e < 1/2 and 

K := < e. (53) 

Then, by Lemma Ul we have that for each i £ Ai there is a state 'kct_-Cm ' 1/ such that 

(we use that log(^) < ^ for 0 < e < 1/2) 


5'max(7rCi'"C'M, IIPCi 


' PCm. ) ^ 


5(r||p) -h2 


(54) 


and 


Then, starting with Eq. ([4^, 


kci-CM, - 111 < 2^/e -h V^. 


1 


E ||tc - Pcili < -pry X] X] 

c&Ci |Ai| , ^ ^ 


\TCi - TTC,-111 


1 


+ _ ^ |£W)|V2 




a,d I A/ 


1 




|Ail . 


E 1 ^. 


hi) 11/2 




< (2 + \/8)Ve 

b,c,j4g,j5^ 



-PC,. II? 



111 


ln(4) 

*S’(vrc,||pc,.) 



S'max(vrCi-C'MjlPl ® ® 


O I ‘I 

<i2 + VS)Vi^ ^4) <'/f 

^ y [n-l + 1 )'^ e 


(55) 


(56) 


^ This is an easy consequence of subadditivity of entropy. Indeed, S{ti-Ai---Am IIpai ® ■ <8) Pa„) = —S{tvai---Am) 

fr(vrAi...A„ log(pAi ® ■■■ ® PA„)) > -E*il S'(7rA^) - tr(7rAi...A„ log(pAi ® ■■■ ® PA„)) = E/1 i ^(TrAjlpA,-)• 
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where we used the (triangle and) Cauchy-Schwarz inequality to obtain the last (first) line. Now, 
for m = 1 we have k = 0 and for m > 2 we have > 1 such that k < e is implied by (see 

Eq. 


2 Ed —zm - (2'^ - ''^n e < -e , 


which also ensures that r > 0 for integer r? Hence, setting 


(57) 


we have that Eq. 
all states r, p. 


r = —I + 




P tzl 

d n e 


implies Eq. (521, which trivially also holds for e > 1/2 as ||re — pc'||] 


(58) 

< 2 for 

□ 


Ernally, using the bound W{z) < ln(z + 1) we find that Eq. ( |T^ implies Eq. (j^, which proves 
Proposition 1^ 


B. Proof of Theorem[T] 


We proof the following stronger version of Theorem]^ 

Theorem 6. Let the canonical state pT (corresponding to a k-local Hamiltonian as in Eq. 0) with en¬ 
ergy density u{T) and specific heat capacity c(T) have (|, z)-exponentially decaying correlations. Let the 
microcanonical state Te^s mean energy such that 

\e-u{T)\ < ^c{T)T^/N (59) 


and energy spread such that 


^ 2df AT\ 

28Afc,^,,,r\/c(r)r2—^ < d < Vc(r)r2. 


Let e > 0 and write 


s = - log 
e 


Vn 


^kg,z,T In^'^(iV) 




lfl£f^,l<l< and s are such that 






s < e(n — I -\-iy 


then 


E lire - pcili < ( V2 + 2 + v^ln(2)) V^. 
CeCi V J 


(60) 


(61) 


(62) 


(63) 


Here, [•] denotes the smallest integer not less than ■ and W the solution to z = VE( 2 ;)e^(^), W(z) >0 
(one of the real branches of the Lambert W function). 


_ Cii 

^ The left hand side of Eq. pTj is lower bounded by (n — Z + l)e «d j such that for n + 1 > 21 we have (I + 
r)err /(^d) > le^/{(.d) > (I — l)err/(^d). 
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This theorem is a direct consequence of Proposition and the following lemma. Theorem 
follows from Proposition the following lemma, and the bound (we use the assumption N > 2, 
the fact that > (c(T)T^)“^/^ = px) — {ir[H and e < 1/2, 

which, as we recall, is w.l.o.g.) 


(64) 

^ In^'^(A^) + (5 + ez) In(jV) 

eln(2) 

Lemma 7. Let the canonical state px (corresponding to a k-local Hamiltonian as in Eq. 0) with energy 
density u(T) and specific heat capacity c{T) have z)-exponentially decaying correlations. Let the state 
T G T>(spfln[{|i/)}i,6Me,J) 




+ ln(iV^+i) 


e-n(T)| < ^Jc(T)T‘^/N 


and 

28Afc,^,,,rV'c(^^^ < <5 < v/c(^. 

Then 


S{t\\px) < -S(t) + log(|Me,5|) + log 


y/N 


(^k,^,z,T In^'^(A) 


,56v^Afe,5,,,Tln2‘*(Ar) 


(65) 

( 66 ) 

(67) 


Proof. We write = NT‘^c{T), p = Nu{T), and define 


Z{T,e,5):= ^ 

Then, for any 5, 5 > 0, any e, e G R, and any state r = u'fsM^ s 
S{t\\px) = -S{t) -tr[Tlog(pT)] 


■S(t)+ r„,„l„g(Am,- 

\MT,e,S) 


Z{T, e, 


< -S’(r) + log 


ZjT) 

-\M~^\ QieN-el^HVN+5VN)/T'\ 


.Z{T,e,6)' y 

We will choose S and e below after Eq. ( [^ and continue with bounding 

Z(T,eJ) 


Z(T) 


^ {v\px\t^) = (i^lPrk) 

v.\E^-eN\<5y/N 

= {^\pt\v) - Y 

i/\ Ei,<eN-\-6y/N v: Ei,<eN—5y/N 

> F{eN + 5 Vn) - F{eN - 5 Vn) 

> G{m + 6 Vn) - G{eN - 6 Vn) - 2 sup \F(x) - 


( 68 ) 


(69) 


( 70 ) 


G(x)|, 














14 


where the (Gaussian) cumulative distribution (G) F is defined in Lemma By the mean value 
theorem, for some x G S\/N), 


G{m + jViV) - G{m - 6 Vn) = 26^/N 




{eN — 
2o^ 


i.e., by Lemma 


Hence, for 


we have 


Z(r,e,J) ^ 2 - ^ 1 (|.iv-.i+^ ln2<^(iV) 


Z[T) a V27r 


- 2A 


2 ~ /— 1 

-6VN—=e 2 ^^ 

u V27r 

Z(T,e,5) 


k,i,z,T- ^ 

Vn 


ln2'^(A) 


Z{T) ^ • 

We now set eN = eN + SVn — 5y/N. Assuming |eA — ^| < cr and 5 < 5 < a/y/N, this choice 
implies |eA — p| + 5y/N < \eN — p| + 6y/N < 2a, i.e., the condition in Eq. ( [7^ is implied by 


(71) 


(72) 


(73) 


(74) 


2 ^ ViV VF 

We return to Eq. ([^ and set 5 = Jq to find 


< 6 < 6 < 


a 


y/N' 


SirlW) < -5(r) + .og . 


(75) 


(76) 


Einally, as we assumed that <5 > So, lE'fc-eAl < Vn5o implies |£^fc-eA| < lE'fc-eAl + leA-eAl < 
\6^/N - 6oVn\ + 6oVN = 5y/N such that |Me, 5 j < \Me^s\. □ 


C. Proof of Corollary]^ 

The first part follows directly from Proposition Lemma and the bound in Eq. 
By Ref. IH^ . for any e > 0, with probability at least 

We,5\e^ 


one has 


1 — 2 exp — 


kc - (a,5)c||i < e + 


ISvr^ 


ds 

<e + 


4 " 


vIBTI 

Eurther, as in the proof of Eemma|^(see the discussion around Eq. ([74|), 

S{qt) = Nu{T)/T + \n{Z{T)) 

( y/N 

< Au(r)/T-eA/T + ln(|Me, 5 |)+ln ' 


Afe,g,,,rln2kA') 

< 2a/T + In (|Me, 5 |) + In (A^) < 2(vk(T) + 1)Vn + In {\Me,5 \), 


(77) 


(78) 


+ 5oVn/T (79) 
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where we used N > 2 and A.k,^,z,T > {c{T)T'^) pt) — {\x[Hpt])‘^) > N 

to obtain the last line. Hence, with probability at least 


we have 


1 — 2 exp 



N s{T) 


2(v^+l) ' 

Vn 


kc - {Te,5)c\\i < e + exp 



2v^ + 2\ 

Vn ) 


(80) 


(81) 


D. Auxiliary Lemmas 


The following auxiliary lemmas were used in the proofs above. The first is the main result of 
IflTl , a Berry-Esseen bound for quantum lattice systems. 

Lemma 8. On A = {l,...,n}^'^ with N = V > 1 sites let H be a k-local Hamiltonian as in Eq. Q and 
let p a state with (^, z)-exponentially decaying correlations. Let 


F{x)= X] h = l^{pH), V = tr {p{H - pf) , 

k: Ef^<x 


and 


G{x) = 


1 


(m-m) 


V^TuH J-oo 


dye 2 < t 2 


the Gaussian cumulative distribution with mean p and variance V. Then 

\V'^{N) 


where 


sup \F{x) — G(x)| < A- 


A = C,(m”<*^-«T+l))'A.ax 


Vn ’ 


1 


1 


(82) 


(83) 


a/VN [max{/c,^}( 2 ; +l)ln(iV)’fj^/A 

and Cd> 1 depends only on the dimension of the lattice. 


(84) 


(85) 


The next lemma was originally proven by Datta and Rermer in IIZTI . in a different formulation, 
and appeared in a form equivalent to the one bellow as Lemma C.5 of [221 . 

Lemma 9. Let it, p, p G be such that S'niax(fr||p) < A and k := 2^\\p — p||i < 1. Then there is a 
state TT such that 


and llvf — 7r||i < V^- 


5'max(7r||/5) < A + log 


1 — K 


( 86 ) 
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Proof. The statement follows from Lemma C.5 of Il22l with Y = 2^p and A = 2^\p — p\. For 
completeness, we give the proof following |22J and I12TI : We have 

^<2^p<2^p + 2^\p-p\=Y + A. (87) 

Let T = (with the inverse the generalized Moore-Penrose pseudomverse) and 

TT = rifrt/tr[rtr7f]. We find TifTt < yV2(y + A)-V2(y + A)(y + A)-V2yi/2 < y_ Further, 
Tty = (y + A)-V2y(y + A)-V2 < i such that 

tr[(i - r'^T)ii-] < tr[(i - r'^r)(y + A)] = tr[y + A] - tr[r+r(y + A)] 

= tr[y + A] - tr[(y + A)-i/ 2 (y + A)(y + A)-i/ 2 (y y 
+ tr[(y + A)-^/2^(y + A)-^/2(y + A)] 

= tr[A(y + A)-^/\Y + A)(y + < tr[A], 

i.e., 

tr[rlTif] > 1 - tr[A] = 1 - K > 0 (89) 

such that 

, = r.rtMTtT.i<y^<^P, (90) 

i.e., 5max(7r||/o) < log(j^). Nowlet 1-0) be a purification of if, tr/i[I (-01] = ff, and write |'0') for the 

unnormalized vector \ f;') = T ® 1|V’)- Then {iplip') = {'flT 0l\ijj) = tr[7fT] such that, as < 1 
(which follows from Tf T < 1), we have 


1 - \{'f\ip')\ < 1 - 




= tr 


if 1 — 


T + yt 


< tr 


(y + A) (1 - 


y+ yt 




= tr [(y + A)] - ^tr ^(y + A)yi/ 2 (y y a )-^/2 
= tr [(y + A)] - tr [(y + A)i/ 2 yi/ 2 l < 


1 


(y + A)(y + A)-^/2yi/2 


(91) 


Finally, 


TT — TT h = 




\'^)W - 


mu’' 

tr[ytyff]_ 


< 


IV')(V’I - 


mm 


tr[ytyif] 


-‘’f- 


(92) 


□ 


Lemma 10. For all pai-Am G T’((C'^)®^) 

M 

Wpa^-Am - Pa^® PAmWi < D‘^^corp^^,„^^{Ai - ■ ■ (93) 

i=2 
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Proof. By Lemma 20 of Ref. 1231 , for every j = 1 ,..., M, 

\\pAr-A^ < D‘^corp^^,„^^{Ai - ■ ■ (94) 

Then by a telescoping sum and triangle inequality, 

M 

IIPAi-Am -PAi ® •••®PAmIIi = \^{L3-Lj-l) ^ 

i=2 

M 

^ J2\\ {pAi-A^ - PAi-Aj.i (8) pAj) ® PAj+i (8) • • • ® PAmIIi 

i=2 

M 

< D‘^^corp^^,„^,{Ai - ■ ■ Aj_i,Aj), (95) 

1=2 

with Lj = pai-a, ® PA,+i ( 8 ) • • • ® pam- 
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